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Abstract 

We present a study of the so called relaxed field equations of general relativity in terms of a decomposition of 
the metric; which is designed to deal with the notion of particles. 

Several known results are generalized to a coordinate free covariant discussion. 
We apply our techniques to the study of a particle up to second order. 
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1 Introduction 

The notion of particle is fundamental to the Newto- 
nian mechanics framework; in fact, this whole theoretical 
framework can be constructed in terms of the notion of 
test particles and massive particles. It is then natural 
to ask whether this notion can be translated to other 
frameworks, as is general relativity. 

Within general relativity one understands Newtonian 
mechanics as the limit of weak field and slow motion. 
So we know that one can regain the notion of particle 
in this regime. Also in general relativity, the concept 
of test particle is a natural one, which allows to discuss 
several physically interesting situations. 

At first sight it is not at all clear that one can extend 
the notion of particles (non-test) to the realm of gen- 
eral relativity. To begin with, if one imagines a process 
in which one shrinks the sizes of an object to obtain a 
point like object, one knows that at some moment in 
the process one would end up with the formation of a 
black hole, which has a characteristic size. However, 
the post-Newtonian approach to compact objects is fre- 
quently constructed in terms of the notion of particles; 
although post-Newtonian systems are normally required 
to have weak fields and slowly moving objects. 

It is interesting to note that the most simple black 
hole, namely the one describing a vacuum spherically 
symmetric spacetime, can be expressed in terms of the 
so called Kerr-Schild decomposition. In this way, the 
Schwarzschild black hole, whose maximal analytic ex- 
tension is described in terms of the well known causal 
conformal diagrams, when expressed in the Kerr-Schild 
decomposition shows a point like description in terms of 
the flat reference metric of the Kerr-Schild form. 



This indicates that it might be possible to give a parti- 
cle notion to a compact object in general relativity when 
expressed with respect to background reference metrics. 

If one intends to study the problem of a systems com- 
posed of several compact objects, it appears as an ap- 
pealing strategy to use approximation techniques for 
solving the field equations. Several problems are related 
to this. 

In building approximation schemes for the study of 
the field equations in general relativity it is often useful 
to recur to the relaxed form of the field equations; that 
we recall below. Also, it frequently useful to decompose 
the physical metric in terms of a background metric. In 
this work we plan to study both techniques. 

In the process of decomposing the metric a key issue is 
the notion of gauge, since in general one has more than 
one way to decompose the physical metric. In order to 
study this issue we bring the techniques used by Friedrich 
in his study of the hyperbolic nature of the gravitational 
field equations. We will present here a generalization of 
Friedrich's results that is convenient for our discussion. 

Although wc work with coordinate independent ex- 
pressions, we also relate our work with the widely used 
harmonic gauge condition; and take the opportunity to 
restate Anderson's result in a coordinate independent 
fashion. 

An approximation scheme is suggested in which the 
previous studies are taking into account. 

We apply our techniques to the problem of a single 
particle up to the second order. 
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2 The decomposition of the met- 



ric 



Let us express the metric g ab of the spacetime M in terms 
of a reference metric r\ a b , such that 



gab — Vab + h a b 



(1) 



Let d a denote the torsion free metric connection of r) a b 
and V a the torsion free metric connection of g ab ; then 
one can express the covariant derivative of an arbitrary 
vector v by 

V a v b = d a v b + T b c v c ; (2) 
and one can prove that 

r a c 6 = \g cd (d a h bd + d b h ad - d d h ab ) = r b c a . (3) 

Let us observe that 
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Fa b9ce — Laefc — ^9° See (d a hbd + 9bh ad — O d h a b) 
= ~ {dahe + b h ae - d e hab) ■ 



(4) 



The relation between L and the curvature tensor can 
be calculated from 



[v a ,v b ]v d = (d a r d c -d b r d c 



+r a d e r b 6 c — v b d e L a e c ) v c 

+ QnbM 



(5) 



abc 
be v i 



R nh d v c 



where O is the curvature of the d a connection. Then the 
Ricci tensor can be calculated from 



Rac — R-abc 

= e oc 

where ac is the Ricci tensor of the connexion d a 



( 6 ) 

p\ I a p o o jit _i_ F " F e f " F e • 



3 Auxiliary functions or gauge 
vector 

Let us consider four independent auxiliary functions x^, 
with fi = 0, 1, 2, 3. Then let us observe that 

,g ab V a V b ^ = 5 ab V a d b ^ = g ab d a d b x» - g ab T a c b d c x^. 

(7) 

Then, if / e denotes the inverse of d c x^, which exists by 
assumption of the independence of the set x^, one has 



.9 Qb r a c b = - ( 3 ab v a v b ^ - 9 ab d a d b x») /; = H^I* 

where we are using 

H* = -g ab V a V b x^ +g ab d a d b x^. 
Alternatively, let us define the gauge vector J^ c 



(8) 
(9) 



(10) 



which implies 

iP = d c x" Jf c = Jf(x>*); 
so that one also has 



g ab r a 



(ii) 



(12) 



These equations show the relation that exist between 
working with a coordinate system, given by the set of 
functions and the gauge vector J4? c ; which does not 
need any reference to coordinate systems at all. In what 
follows we will try to use the covariant approach that 
employs the use of the gauge vector Jj? c . We empha- 
size that Latin indices are abstract; and therefore our 
expressions are coordinate independent and covariant. 

Then, the Ricci tensor can be expressed by 
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Rac - @ac + (©fea/^ec + ^bed^ea + 2<d bca e h ed ) 

+ \g bd d b d d h ac - d (a (g c)e J? e ) + g ed T a e c Jf d 



9 bf 9ed^af^bc ~ 9 ( r a Mr &cd + F bd T bad ) 



(13) 



Let us note that if the vector field M" 1 is given by ( 12 ), 
then for any function x^ one has 



g ab V a V b x» + g ab d a d b x» = Jtf(x»). 



(14) 



In the standard studies on approximations to the solu- 
tion of the field equations, one frequently finds the choice 
of harmonic coordinates for the set of the x M 's; however, 



in equation ( 13 ) one can see that only the vector field 



Jf? c appears, without any reference to a choice of aux- 
iliary functions. Therefore one could just refer to the 
gauge vector Jt? c . 

4 The field equations in relaxed 
covariant form 

Previous to the discussion of the relaxed covariant 
field equations, we would like to refer to the work of 
Friedrich[3] and its extension to this coordinate inde- 
pendent discussion. 

4.1 Priedrich's theorem without the use 
of coordinates 

The field equations are 

II, Stt/s ( T ac - l -g ac g bd T bd J . (15) 



Equation (3.22) in reference [3] can be obtained from 



( 13 ) by expressing it in a coordinate frame and neglecting 



the terms. In this way, one would obtain the analogous 
expression where all the appearance of d derivatives are 
replaced by coordinate derivatives c^, the tensors F are 
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replaced by the Chrirstoffel symbols and one uses F e = 
—H e ; namely: 
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g p dpd s g aa +5 £ ( CT V Q) F £ 
f*g e5 T%Tl a - \ I> aJ + T^T PaS ) (16) 



-87TK ( T aa - -g aa gP 5 Tps 



Friedrich has studied [3J this system introducing the 
notion of "coordinate gauge source" 



(17) 



Subsequently, Friedrich studied the case in which is 
given arbitrarily. 

Then, we can rephrase Friedrich's theorem in the fol- 
lowing form: 

Theorem 4.1 Let x^ be four independent functions that 



are used as a coordinate system. If is a solution of (16) 
together with the matter equations such that on the initial 
surface one has = W, VpF* = V^VV^x", 
then g^ v is in fact a solution of Einstein's field equations. 

This theorem can be understood in two ways. In one of 
them, we think that the four coordinates x^ are given 
and then the theorem checks whether the F M 's satisfy 
the above equations. In the other way, one think that 
the F^'s are given and then the theorem checks whether 
there exists a coordinate system of x^'s such that the 
equations in the theorem are satisfied. 

From the fact that = — F M , one deduces, using the 
same techniques as in [3J, that the generalized Friedrich's 
theorem holds, namely, consider the four functions 
as given a priori, then: 



Theorem 4.2 If g a b is a solution of ([15]), with the de- 
composition of the metric as in ([T]) and with the Ricci 



tensor as given by (13) with J4? c = H^I^, together 



with the matter equations such that on the initial surface 

one has = -(g ab V a V b x» + g ab d a d b x^), V C H» = 
-V c {g ab V aV bx* 1 - g ab d a d b x^), where are four inde- 
pendent scalars, then g ab is in fact a solution of Einstein's 
field equations. 



This result gives great freedom in the problem of find- 
ing solutions of the field equations in terms of a refer- 
ence metric. Suppose that one solves (15 1 for a given 
vector field H^. Also assume that one can solve for the 
functions x^ such that g ab V a V b x^ = -H*. Then, let 
us build a flat metric 77 so that g ab d a d b x^ = 0; which 
in particular can be satisfied if the cc^'s are thought as 
Cartesian coordinates of 77. In this way one would ob- 
tain H^I^ = g ab T^ b , and so have a solution of the field 
equations. 

It also might be of interest to researchers in numeri- 
cal relativity, since it provides the possibility to use any 
coordinate system; i.e., not necessarily an harmonic one. 

Instead, one could have a proposition that does not 
refer to the auxiliary functions whatsoever; namely 



Theorem 4.3 If g ab is a solution of (15), with the de- 
composition of the metric as in and with the Ricci ten- 
sor as given by (13), together with the matter equations 



such that on the initial surface one has = </ r Q c b , 
V ' d-^ c — ^d{g ab ^ab)' then g ab is in fact a solution of 
Einstein's field equations. 

This theorem can be understood in two ways. In one 
of them, we think that the metric rj ab is given and then 
the theorem checks whether the vector Jif c satisfies the 
above equations. In the other way, one think that the 
J(? c is given and then the theorem checks whether there 
exists a metric rj a b such that the equations in the theorem 
are satisfied. 

4.2 Relaxed covariant form of the field 
equations and a generalization of 
Friedrich's theorem 

Alternatively one can use the form of the field equations 
in terms of a slight different logic. 

When we use the expression of the Ricci tensor as 



given by (13) in (15), without assuming that .34?° is 
g abT ab> namely 
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g bd d b d d h ac - d {a {g c)e ,^ e ) + g ed T a d c Jf e + 9 fl 



- g bf ged^af^bc - 2 (rj^r&cd + r c bd r ba(i ) 



8itk [ T ac - -g ac g oa T bd ) ; 



(18) 



we will refer to these as the relaxed field equations [4] , 

Using the standard harmonic gauge technique, one 
would say: solve the relaxed field equation in the coordi- 
nate frame, with = 0, and then require the equation 



g bd V b V d x» = 0. 



(19) 



In the standard approach one makes use of coordinate 
basis; therefore the previous statement would be the 
complete story. However in our case, H 11 has a second 
term where two covariant derivatives of x^ with respect 
to the metric 77 appears. At this point it is important to 
notice that if we have the solutions x^ 1 from ( 19 ) then, 



on constructing r\ with this as a Cartesian coordinate 
system, one would obtain ii 7 ' = 0. 

In some occasions it is preferable to work with a dif- 
ferent set of equations. In this respect, several authors 



have indicated that actually to request equation (19) is 
equivalent [213] to demand 



0. 



(20) 



When dealing with Einstein equations in the relaxed 



form, and treating the vacuum case, equation (20) is 



understood as the condition that the divergence of the 
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Einstein tensor must be zero (which of course is identi- 
cally zero in the non relaxed form). 

Let us study the relation between the divergence of 
the energy-momentum tensor and the vector J^ e . One 



can write the relaxed field equations ( 18 ) in the usual 



form in which on the right hand side we have just the 
standard term — 87rreT ac ; and so on the left hand side, 
the terms involving M' e would be 

- V (a (g c)e ^ e ) + \gacg ef V e .^ f - (21) 

where we have used that the term 8ttkT contributes with 
the term 

-. 9 e/ V e ^. (22) 

Then in taking its divergence, on the left hand side, 
the terms involving J^f e are 

3 ah V a (-V (o (g c)e J? e ) + \g a cg ef ye^ ■ (23) 

If we replace Jf e by T e = g ab T^ b , the divergence of the 
left hand side would be identically zero, since the Ein- 
stein tensor has divergence zero. Therefore we conclude 
that the divergence of the stress energy-momentum ten- 
sor is 



8n K g ab V a T bc = 5 ab V a ( -V (6 (g c)e (J? e - T e )) 
+ lg bc g de V d (g e f(Jf f -T*)) 



(24) 



Therefore, the stress energy-momentum is conserved if 
and only if 



(25) 



g ab Va[ -V (b ( 5c)e (^ e -r e )) 

+ \g b cg de V d (g ef (Jf f - r/))) = 0. 



Working out the relations, one finds that the previous 
equation can be expressed as 



o = -^ ce g ah v Q v h (j? e - r e ) + l -R ce {^ e - r e ). 



(26) 



Which coincides with Friedrich calculation. 



It follows that if one solves equation ( 26 ) such that on 
an initial hypersurface 34 ob = T b and V a J4? b = V a r b , 
then the energy-momentum tensor will be conserved in 
the evolution of the system. Furthermore, one also de- 
duces that: 



This is a rephrasing of Friedrich's theorem applied to a 
decomposition of the metric and to its general relaxed 
covariant form of the field equations. 

It is interesting to remark that Andersonjl], using a 
retarded integral expression for h, was able to prove 
the equivalence between the conservation of the energy- 
momentum tensor with the harmonic gauge condition. 
In relation to this let us remark that if the set of func- 



tions x^ is obtained from the solutions of (19); and one 



uses them as harmonic coordinates of the metric 77, then 
one deduces that r c = 0. And also, if T c = 0, then 
Cartesian coordinates of 77 are harmonic coordinates of 
g. This means that we can state Anderson's result in a 
coordinate independent way, namely: 

Theorem 4.5 Let h be the retarded solution, with respect 
to a flat metric 7/ of the relaxed field equations together with 
the matter equations of state, such that Jf? b = 0, then the 
conservation of the energy-momentum tensor implies that 
g a t is a solution of Einstein equations. 

5 The approximation method 
and the treatment of particles 

The approximation method that we introduce below, is 
adapted to the treatment of particles; therefore, it is 
convenient to begin by treating the problem of one single 
particle in the context of linearized gravity, in order to 
clarify some of the techniques. 

5.1 The gravitational field from one par- 
ticle in linearized gravity 

5.1.1 The description of a particle 

Let us consider a massive point particle with mass rriA 
describing, in a flat space-time (M ,?^), a curve C 
which in a Cartesian coordinate system x a reads 



z»(t), 



(27) 



with r meaning the proper time of the particle along C. 

The unit tangent vector to C, with respect to the flat 
background metric is 



(28) 



that is, r](u,u) — 1. Now, for each point Q — z(t) of 
C, we draw a future null cone Cq with vertex in Q. If 
we call Xp the Minkowskian coordinates of an arbitrary 
point on the cone Cq, then we can define the retarded 
radial distance on the null cone by 



Theorem 4.4 If, given the metric 77, one solves the 
relaxed field equations for h together with the matter 
equations, which include the conservation of the energy- 
momentum tensor, such that Jff b = T b and V a Jf b = 
V a r b on an initial hypersurface, then g^ is a solution of 
Einstein equations. 



(29) 



The energy momentum tensor of a point particle 
is proportional to mv a Vb] where m is the mass and v a 
its four velocity. We are distinguishing between the 
unit tangent vector u a and the four velocity vector v a , 
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because in future works we would like to consider the 
possibility to normalize the vector v with respect to a 
different metric. In order to represent a point particle 
must also be proportional to a three dimensional 
delta function that has support on the world line of the 
particle. 

We will suppose that the particle does not have multi- 
polar structure. Then, given an arbitrary Minkowskian 
frame (x , x , x 2 , x 3 ), we will express the energy momen- 
tum by 

T(°)° b (x =z (T ),x\x 2 ,x 3 )=m A v a (T )v b (T ) 

8(x 1 - z^Tp^Sjx 2 - z 2 {t ))S{x 3 - z 3 (r )) (30) 

5.1.2 The first order solution 

The retarded solution, in terms of Green functions, for 



the relaxed field equations ( 18 ) for particle A, in which 
we take J(f b = and r\ a flat metric, is 



so that in general 



7 ( ? 

Jab 



= 1 



-Am a 



2niA 



V a Vb 



Vab 



~V a Vb- 



(31) 



(32) 



In these equations we have considered the definition 

Va = VabV b . (33) 

It is interesting to realize that the exact inverse of this 
metric is 



(l)ab 



1 



4m £ 



l 



1 



(34) 



Note that one can solve for h ab for an arbitrary motion 
of the particle; however, the complete solution of the 
problem involves having to set also T b = 0; which in 
terms of a coordinate frame treatment is equivalent to 
the harmonic condition. Then, recalling, as mentioned 
previously, that Anderson has proved [T] the equivalence 
between the harmonic condition and the divergence free 
condition on the energy-momentum tensor; one deduces 
from this, that for the case of the energy-momentum 
tensor of a particle it implies its geodesic motion. 

5.2 Iterative approximation method 

Now we present a general iterative method to solve the 
relaxed field equations. 

First of all, let us note that given the decomposition 
lb and defining the tensor h ab from 



where 



(35) 
(36) 



that is rj ab is the inverse of rj ab , one can always express 
the inverse g ab in the form 



g ab = 7 fb _ ~ h ab _ d ab_ 

Then making the contraction 

g ab g bc = S a c - h ab h bc - d ab g bc = 6 a c ; (38) 

one finds 



(37) 



d ab = -h ad h dc g cb - 



(39) 



which can be considered an implicit equation for d ab ; but 
it also shows explicitly that d is quadratic in terms of h. 

This suggests the natural series d%, da, d^, d$, d^, ... de- 
fined by 

df = -h ad h dc ( V cb ) , (40) 



df = -h ad h dc (r] cb -h cb 



iab uad l. ( ^cb 7cb jcb 

d n = -h h dc 1 77 -h - d (n _ 2) 



(41) 



(42) 



for natural numbers n > 3. It is clear that d n is order 

However, we have seen that in the first order solution 
for a single particle, the inverse of the metric has a term 
which is conformal to the flat metric 77; which it will 
be convenient to take into account. For this reason we 
propose the following method of approximation where 
this issue is considered. 



The idea is to express ( 18 ) and eventually ( 19 ) in the 
form 

<prf a, d a d b f = S(f); (43) 

where <p rj ab is the term proportional to r] ab that is con- 
tained in g ab ; while the general case would be to consider 
just rj ab d a dbf for the left hand side. This equation can 
also be expressed by 



V ab d a d b {v f) = 8{<p f) + S(f); 



where 



8(<p f) = V ab dad b & f)- V ^"dadbf. 



(44) 



(45) 



Now one would like to solve equation ( 44 ) by iterations 



Let us define the sets such that for j = 0, one 
takes h — 0, x M to be harmonic functions of the metric 
r\ and ip = 1; and for j > 0, / < - J - > is the solution of 

rfdadbipU-V /«) = s(<pU-V /O-D) + 

(46) 

using the retarded Green function. As we have seen 
above, ip^ 1 ' clearly arises in the first order calculation. 

The application of this method to the first order, for 
a single particle, reproduces the calculation explained 



in subsection |5.1.2| Next we study this case at second 
order. 
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5.3 The second order solution 

Let us remark that the first order solution is stationary 
and spherically symmetric. This structure transports to 
the second order solution. 

(2) 

The equation for h y ab ' is 



9 9 e d l a f 1 bc 

1 /Wl)6d r (l) i r (l)hd r (l) 
2\ a bcd^~ 1 c L bad 

87TTI 1 ). 



(47) 



We will call the right hand side, the tensor T; which 
has the structure 

T = a{r)dt 2 - P(r) (dxf - 7 (r) (x ■ dxf ; (48) 

where we are using the three dimentional notation x = 
(x, y, z) and where 



2m 2 



o 



(2m + r) 2 (2m — r)r 

— 2m 2 (4m + r) 
(2m + r) 2 (2m — r)r 2 ' 
8m 2 (2m 2 + mr — 2r 2 ) 



inmS 3 (x), (49) 
(50) 
(51) 



' (2m + r) 2 (2m- r) 2 r 4 ' 

Therefore one assumes for the same form, namely 
h {2) = A(r)dt 2 ~ B(r) (dxf - C(r) (x ■ dxf . (52) 



In this way one has 



,(2) _ 

,(2) _ 

Hi ~ 

,( 2 ) - 



.4. 



-(B + (x*) 2 C), 
-x l x 3 C\ 



(53) 
(54) 
(55) 



where the index i,j — 1, 2, 3 denote spatial coordinates. 
One can see then that the equations to solve are 



□ C- 



□ A = 
4dC 
r dr 

□ B = 2(13 + 0); 



2a, 
= 27, 



(56) 
(57) 
(58) 



where we are using the symbol □ to denote r] bd dbdd- 

We can solve these equations in two ways, either using 
Green function techniques, or, recalling the stationary 
nature of the solution, just integrating the Laplace oper- 
ator. For this presentation we choose the second option. 
Let us note that for any function f(r) one has that 



Uf(r) = -V 2 f = - 



1 d 

r 2 dr 



dr 



(59) 



Therefore one can find A by two consecutive integra- 
tions, obtaining 



1 



A(r) = A (l- — )ln 



2 m 



ka 2 - 



kai 



(60) 



Similarly one can see that the function C(r) satisfies 



l_d^ f 6 dC . _ ^ 
r 6 dr \ dr 



which after integration gives 

1 / 9m 

°( r ) = 7n 2^5 ~( r + 288 ™ 5 ) ) 

40m z r° \ r 

9m 

+ (416m 5 -7r 5 )ln(l + ) 

r 

+ 128m 5 ln(r) 

- 8fccim 2 + 40fcc 2 m 2 r 5 - 352m 4 r 
+ 16m 3 r 2 - 16mV + 12mr 4 

Then the function B(r) is given by 
1 



(61) 



(62) 



B(r) 



120m 2 r 3 

(288m 5 - 20m 3 r 2 - 30mV + r 5 ) 
, , 2m N 

Hi- — ) 

+ (-416m 5 + 20m 3 r 2 - 130m 2 r 3 + 7r 5 ) 
2m 

ln(l + ) - 128m 5 ln(r) 

r 

- 120fcfeimV + 120kb 2 m 2 r 3 + 8k Cl m 2 

- 40kc 2 m 2 r 5 + 352m 4 r 

-1024m 3 r 2 - 12mr 4 



(63) 



Our choice for the integration constants is: ka\ 



ka 2 — 0, kb\ = — f5 m : kb 2 = 



17 

15 > 



kci 



74 „ 



= m, 
and 



kc 2 = 0. This choice is made taking into consideration 
the exact solution and the integration of the solution 
coming from a Green function approach; that we will 
not discuss here due to considerations of space. By the 
exact functions we mean the metric components of the 
Schwarzschild spacetime in harmonic coordinates; which 
are 



A-exact 
B * exact 
O exact 



1 - 

TT 



= i 



m 
r 



i m 



(64) 
(65) 
(66) 



A graphical comparison with the exact functions of 
the Schwarzschild solution in harmonic coordinates are 



shown in figure 5.3 



It can be observed that in second order one obtains 
an excellent comparison of the solution with the exact 
values of the metric components; even for very small val- 
ues of the radial coordinate. Although this comparison 
has limited value, it is in any case remarkable that it is 
only necessary to go only to second order to obtain such 
a good approximation. 
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Figure 1: Comparison of the functions calculated in the 
second iteration with the exact values. 

6 Final comments 

We have presented an study of an approach to the grav- 
itational held equation through the relaxed covariant 
form of them. The whole approach is intended to deal 
with the notion of compact objects. 

The relaxed field equations was studied using Friedrich 
approach to the problem and we have also refer to An- 
derson's result in the field of harmonic conditions. 

We have generalized Friedrich results to a covariant 
formulation in terms of a decomposition of the metric. 

Anderson's result has been restated in a form that 
does not make reference to coordinate conditions. 

We have presented an approximation method that can 
be applied to the notion of particles in general relativity; 
and which is successful in second order for the case of a 
solitary compact body. 

It is our intention to apply these techniques to the 
problem of a binary system in general relativity. 
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